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^ ; Abstract 

The Nambu- Jona-Lasinio model (NJL) is extended to incorporate heavy 
. quark spin-symmetry. In this model baryons containing one heavy quark 

are analyzed as bound-states of light baryons, represented as chiral solitons, 
and mesons containing one heavy quark. From related studies in Skyrme 
type models, the ground-state heavy baryon is known to arise for the heavy 
meson in a P-wave configuration. In the limit of an infinitely large quark 
mass the heavy meson wave-function is sharply peaked at the center of the 
chiral soliton. Therefore the bound state equation reduces to an eigenvalue 
problem for the coefficients of the operators contained in the most general P- 
wave ansatz for the heavy meson. Within the NJL model a novel feature arises 
from the coupling of the heavy meson to the various light quark states. In this 
respect conceptual differences to Skyrme model calculations are discovered: 
The strongest bound state is given by a heavy meson configuration which is 
completely decoupled from the grand spin zero channel of the light quarks. 
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I. INTRODUCTION AND MOTIVATION 



Recently, there has been considerable attention given to properties of hadrons containing 
a single heavy quark with mass, Mq, much larger than the typical scale of strong interactions, 
Aqcd- Henceforth these particles will be referred to as "heavy hadrons". These "heavy 
hadrons" are frequently described within the heavy quark effective theory (HQET) 
which represents a 1/Mq expansion of the heavy quark content of QCD. The heavy quark 
effective Lagrangian is set up from the heavy quark transformation 

Q { } s \x) = ^ exp (iM Q v ■ x)Q(x), (1) 

which disentangles the large (I) and small (s) components of the heavy quark spinor, Q(x). 
In eq ([]]), refers to the velocity of the heavy quark spinor which, e.g. in its rest frame 
takes the form, t> M = (1,0). The Dirac operator for the large component is expanded as 
v ■ D + 0(1/ Mq). Since the leading order term commutes with the quark spin operator, 
heavy mesons differing only by their spins, as e.g. pseudoscalar and vector mesons, are 
degenerate in the heavy quark limit, Mq — ► 00. The empirical spectrum of heavy mesons 
@, m D = 1.87GeV « m D * = 2.01GeV and m B = 5.28GeV « m B * = 5.32GeV represents 
the most convincing evidence of the heavy quark symmetry^. This degeneracy is in contrast 
to the spectrum of mesons containing only light quarks. Here the pseudoscalars are much 
lighter than the vector mesons. This is understood by interpreting the light pseudoscalar 
mesons as the would-be Goldstone bosons of spontaneous chiral symmetry breaking. 

The incorporation of the heavy quark symmetry into effective meson theories is accom- 
plished by introducing the heavy meson field, 7if h , in the heavy quark rest frame Here, 
I and h denote the light and heavy quark spins respectively, while a = u,d refers to the 
isospin of the light quark. This field contains pseudoscalar (P a ) and vector (V®) mesons on 
equal footing [Q, 

HI = iP a 5 lh + a\ h V?. (2) 
According to the heavy quark symmetry the heavy spin index (h) decouples completely in 



See for further manifestations of the heavy quark symmetry. 
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effective theories^. The interactions of Tlf h with the light mesons proceeds via the light quark 
indices (I, a) and is governed by chiral symmetry. Such effective models have been employed 
to not only study the properties of heavy mesons, (spectrum, form factors, weak decays) 
but also the spectrum of heavy baryons [|5|-|T0|j. In this context the soliton picture for 
baryons is commonly adopted. Within this approach the baryon number is carried by a 
solitonic configuration of the light meson fields. Heavy baryons emerge when a heavy meson 
is bound in the background of this soliton. This picture relies on the bound state approach 
to strange baryons initiated by Callan and Klebanov [Tl||l^]. The main difference here is 
that for the heavy meson both the pseudoscalar and vector meson degrees of freedom have 
to be taken into account as a consequence of the above mentioned degeneracy. A continuous 
transition from the hyperons to the heavy baryons containing c- or 6-quarks instead of an 
s-quark suggests a P-wave configuration for the bound meson. 

In the heavy quark limit the heavy meson wave-function, 7if h , is confined to the center 
(r = 0) of the soliton. For a quantitative determination of this effect it is necessary to take 
into account the interaction amongst the heavy quark degrees of freedom before performing 
the transformation ([[]) ||. However, to leading order in the 1/Mq expansion it is consistent 
to investigate the interaction of a strongly peaked heavy meson wave-function and the light 
quark degrees of freedom while henceforth discarding that part of the action which only 
describes heavy quark degrees of freedom. Up to 0(1/ Mq) corrections the binding energy 
of the heavy baryon is thus determined by the properties of the soliton at r « []§~"0] . In 
that case the bound state equation reduces to an eigenvalue problem for the coefficients of 
the operators in the general P-wave ansatz for 7if h . 



Here we study such bound systems within the Nambu-Jona-Lasinio (NJL) model [13 



extended to incorporate heavy quark symmetry ||14||P|. In contrast to Skyrme type models the 
soliton in the NJL model contains components of non-vanishing grand spin. The grand spin 
is the vector sum of total angular momentum and isospin, G = J + r/2. The appearance of 
additional modes is due to the quark sub-structure of the soliton and leads to a novel coupling 
scheme between the heavy meson and the soliton as will become clear in the course of this 



2 Note that in the rest frame of the heavy quark the time component of the vector field is suppressed 
in the heavy quark limit by virtue of the Proca-type equation of motion. 
3 Similar studies have been carried out in ref [15]. 
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paper. In particular, the largest binding energy is obtained for a heavy meson configuration 
which completely decouples from the grand spin zero channel of the light quarks. 

In section 2 we review the extension of the NJL model to include heavy quarks. In 
section 3 the bound state wave-function of the heavy meson in the background of the 
soliton in the NJL model is constructed. The numerical results are presented in section 4 
and concluding remarks are contained in section 5. Some technical details are relegated to 
an appendix. 



II. HEAVY QUARKS WITHIN THE NJL MODEL 

In NJL-type quark models the interaction Lagrangian is given by a four quark contact 
term. The heavy quark symmetry is incorporated by only taking into account the large 
component Q v =Q$ which is defined in eq. (|1]). Furthermore, that part of the Lagrangian, 
which only contains heavy quarks, Q, is discarded. We refer to ref |TJJ for the explicit form of 



the four quark interaction between the heavy and light quarks. Following the path integral 



bosonization procedure JRl, a composite field, H, is introduced such that the Lagrangian is 



bilinear in the quark fields, 

C hl = Q v (ivd)Q v - Q v Hq - ~qHQ v + -i- tr(M), (3) 

/Ct3 

where q denotes the Dirac spinor of the light quarks in the chirally rotated representation 
q = Qq ||17|| . Denoting by £ the root of the chiral field, U = in unitary gauge, we have Q = 
£ + + k — 6^) 7s- The second and third terms in eq. (§) represent the interaction between 
the heavy and light quarks which is mediated by the heavy meson field, H = (H u , H d ) Q 

H a = hl + m^P a + f a ), H a = loH^ ^ 0) ( ° V a = u,d. (4) 

2 \n a oj 

Finally the last term in eq. (0) is purely mesonic. The coupling constant G% stems from the 
four quark interactions between the heavy and light quark fields. In eq. (|3]) only the leading 
terms of the 1/Mq expansion have been maintained. Therefore, the kernel (iv-d) for Q v is 
that of a massless non-relativistic particle. Integrating out the fermion fields yields 

A = A T + A l m + A h m) 
At = Tr A log D, 

A h m = / d±x tr (HH) . (5) 



The functional trace stems from the path-integral over all quark fields and involves the 
inverse propagator, 



D 



V 



-H iv ■ d 



(6) 



which acts on the light-heavy spinor 




. Here m denotes the light quark constituent 

mass. As a consequence of the chiral rotation and the omission of the light (axial) vector 
mesons, the light meson fields only appear in A l m . Its explicit form may e.g. be found in 



ref. [17]. Finally, p, refers to the Dirac operator for the light quarks, u,d and describes 
the interaction between these quarks and soliton constructed from the light mesons which 
in turn are quark composites. The explicit form of p l will be presented in section 3. 

In eq. (|5|) we have already indicated the need for regularization. For the ongoing dis- 
cussion we will employ Schwinger's proper-time regularization |f[Sjl . In the NJL model this 
implies a continuation to Euclidean space (r = rr 4 = ixq). This defines the Euclidean Dirac 
operator JD E and the Euclidean action A)p = TrAlogD E . Subsequently, its real part is 
replaced by a parameter integral 

r°° ds 



1 + I f°° ds / 
A* = -Tr A logD E Dk = -- / — Tr exp -sD^ 

2 2 Jl/A 2 S v 
with the high momentum contribution chopped off. The imaginary part, 



Ajp = -Tr log Dp CD 



t 



(7) 



(8) 



is finite and remains unchanged. After computing the functional trace the action is continued 
back to Minkowski space. To adjust the model parameters we demand the physical pion 
mass and decay constant, m n = 135MeV and f n = 93MeV, respectively. Then all but 



one parameter of the light quark sector are determined [ 16.19] . For transparency this is 
commonly chosen to be the constituent quark mass, m. Subsequently, the heavy quark 
coupling constant, G3, is fitted to the 5-meson decay constant fs ~ 180MeV as estimated 
from lattice QCD |F2D[ and QCD sum rules |[21|| . Since the present treatment closely follows 
the computation of ref. 0] we dispense with further details although one remark is in order. 
As mentioned above, the light (axial) vector mesons are omitted in the present study. As a 
consequence of the missing n — a± mixing the relation between f n and A is modified whence 
reducing the numerical value of A. This also effects the prediction for G3. We display the 
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resulting values in table |. For the system without the light (axial) vector mesons we have 
re-evaluated the binding energy, Em = m — AM, of the heavy meson by determining the 
root, AM of the Bethe-Salpeter equation for the heavy meson field H. For simplicity we 
expand this Bethe-Salpeter equation in powers of v ■ p, with p^ = — MqV^ being the 
residual four-momentum of H . The coefficients in this expansion depend on the proper-time 
cut-off, A. It turns out that the expansion in v ■ p converges less rapidly than in the model 
with light (axial) vector mesons. An expansion up to at least cubic order proved necessary 
to gain reliable results. In table | we present the predicted root, AM, and the binding 
energy, Em- These results were obtained from an expansion up to fifth order in v ■ p. The 
analytic form of this expansion is provide in the appendix cf. eqs ( p4f) and fl25l) . Apparently 
the binding energy increases with the constituent quark mass. It should be stressed that 
these results correspond to the leading order of the 1/Mq expansion. 

Eq. d6|) shows that only the residual momentum, p M , is affected by the regularization 
as a consequence of the heavy quark transformation ([!]). This is physically meaningful and 
crucial for a consistent interpretation of the model. Indeed, the residual momentum of the 
heavy quark in a hadron containing a single heavy quark arises entirely from its interaction 
with the light degrees of freedom (in this case light quarks) and is thus cut off at the same 
scale A. Note that the cut-off A is significantly smaller than the c- or b- quark masses (cf. 
table |). 



III. BARYONS WITH A HEAVY QUARK 

Before describing heavy baryons as bound systems it is appropriate to briefly explain 
the emergence of the soliton. When integrating out the light quarks an energy functional 
for the light meson fields can be extracted |2^|T^] by first defining the one-particle Dirac 



Hamiltonian, h, from if3p l — idf — h , cf. eq @. Adopting the hedgehog ansatz for the 
chiral field U = exp (it ■ fO(r)) this Hamiltonian reads h = cx-p+m(3 (cosG + 275 r • f sin0). 
The total energy functional is expressed in terms of the eigenvalues, e M , of h. Formally it is 
composed, respectively, of valence, vacuum, and meson contributions 

E tot [Q] = Ncri^M + f°, *- Stl + I d 3 r (1 - cos 0(r)) . (9) 

2 ^ Jl/A 2 V47TS 3 J 

Here Nq = 3 is the number of color degrees of freedom. Furthermore, r] WSb \ denotes the 
occupation number of the valence quark level. This level is defined to be the one with 
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the smallest module le^. Its occupation number is adjusted to describe a configuration 
possessing unit baryon number, i.e. r] vau \ = 1 + (1/2) sgn(e M ). Finally the chiral angle, 
0(r) is self-consistently determined by minimizing E to t p3| , |i~9f . 

In order to compute the binding energy of the heavy meson in the soliton background 
the bosonized action (^) is first expanded up to the quadratic order in the meson field, 
H(r,t) = J(du/2ir) H(r,u) exp(-iwt) 

A% = /M^ tr (M) - J^fl ™ M -U^PHMMR^A (10) 
J 2ir { iV c G3 v ' to + e va i ~ J 

where the arguments (r, u) of the heavy fields have been omitted. In analogy to the static 
energy functional @ this part of the action is composed of a purely mesonic piece as well 
as valence and vacuum contributions. The regularization function, R\ is obtained as an 
expansion in the frequency to 



Ra{^, e) — — ( 1 — sgn(e) erfc 



2e 



— zr^il — sgn(e 



+— \\\- sgn(e) erfc 



2e 2 V 

6 

A 



2A 



e-V* +•••■ (ID 



This expansion is equivalent to the one in terms of v-p for the meson sector. For completeness 



we list Ra up to fifth in the appendix, cf. eq (|27|) . It is not surprising that this result for 
can equivalently be obtained from the expression found for the fluctuations of kaon 



fields in the background of the NJL model soliton |M| once the eigenenergies for the strange 
quark levels are set to zero. We note that in the unregularized formulation only the negative 
energy states contribute to the vacuum part of because 

Ita^^ii^fe). (12) 

In the many body interpretation of the functional integral, A^, one would phrase this result 
such that in the limit A — > oo only the occupied quark orbits couple to the heavy meson. 
The unregularized expression ( |T2"D may independently be derived when considering 

AgWeg.) = Tr {(ip ^ HG V H} (13) 

in Minkowski space. Here 

G v (t,t>) = (t\ (id,)- 1 \t>) = lira / ^ e — — (14) 
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denotes the non-relativistic propagator in the heavy quark rest frame. However, the reg- 
ularization is unavoidable since J2a Roo(^, £/x) diverges logarithmically. Nevertheless, the 
limit fll2|) serves as an independent check on our computations. From the expression ( |i2|) 
it also becomes apparent that the total action fllQ|) is continuous as the energy eigenvalue 
associated with the valence quark orbit changes its sign. 

Next we study the coupling of a grand spin |, P-wave heavy meson to the soliton. Such 
a configuration can be expressed in terms of three radial functions |§ 

H? h (r, u) = r- r ab {u A (r, u)5 lh 5 bc + u B (r, u)r bc ■ a lh + u c {r, u)r ■ r bc r ■ a lh ) X c . (15) 

Here \ a denotes a (constant) isospinor. One might have expected the appearance of three 
different isospinors. However, such an ansatz is waived because it is not an eigenstate of the 
grand spin projection operator. Skyrme model studies || also indicate that the light grand 
spin, which is defined as the grand spin of the operator multiplying x, has to vanish for the 
meson configuration with the largest binding energy. In the heavy quark limit the binding 
energy of this heavy meson is determined by the potential it experiences at the center of 
the soliton, r = 0. In order to extract this potential it is convenient (and consistent) to 
adopt radial functions in eq ([15]) which are non- vanishing only at r = 0. Consequently, the 
heavy meson field only couples to quark levels, which have non-vanishing radial functions at 
r = 0. As can be seen from the decomposition (^) and the coupling scheme exhibited in the 
Lagrangian (|^), these non-vanishing radial functions must occur in the lower components 
of the light quark fields. At first sight this appears to be surprising because the lower 
component of the valence quark eigenfunction happens to vanish at r = 0. However, one 
has to be aware of the fact that the heavy meson couples to the light quark states via 
the chiral rotation, which is off-diagonal at the origin in the presence of the soliton since 
9(r = 0) = -7T, 

(0 it ■ r\ 
(16) 
it ■ r / 

allowing the valence quark to couple to the heavy meson. In addition to the channel with 
grand spin zero, which includes the valence quark orbit, only the channel carrying grand 
spin one contains states, which do not vanish at r = 0. 

Carrying out the angular integrals as well as the flavor and Dirac traces finally yields 
the action as a function of the values which the meson wave-functions in eq (|I5|) assume at 
r = 0. Denoting these by A(u), B(uj) and C(u), respectively, we find 



S 



A 



(2) 
Ih 



duj 
2n~ 



[x ] x) (A*(u),B*{w),C*{u)) Q(u) 
The inverse propagator is a 3 x 3 matrix 



/ A(u) \ 

B(u) 
\C{u)j 



(17) 



in 



N C G 3 



(1 





o\ 




/ 1 -3 -1\ 




/3 


3 


1\ 





3 


1 


+ /oH 


-3 9 3 


+ /iH 


3 


3 


1 


U 


1 






[-1 3 1 J 




u 


1 


37 



(18) 



Denoting by g^(r) the radial parts of the upper components of the light quark wave- 
functions, which are eigenstates of h, the spectral functions are expressed asQ 

|2 



n Jjt^)L + £ ^(o)! 2 ^^) 



At,G=0 



^ |^(0) | 2 i? A (w,e M ) 



(19) 



The sums refer to distinct quark grand spin (G) channels. 

The equation of motion for the heavy meson reduces to a homogeneous matrix equation 
for A, B and C . Its solution requires adjusting u to u such that Q{ojq) possesses vanishing 
eigenvalues. It is easy to show that there exist two distinct solutions to this problem 

2tt 



(1) : 

(2) : 



N C G 3 
2w 

N C G 3 



+ f Q (u ) = 0, with (A, B, C) cx (1, -1, 0) , 
+ h (wo) = 0, with (A, B, C) oc (3 + C, 1, C) 



(20) 
(21) 



The vectors of set (2) obviously span a two dimensional vector space causing the associated 
eigenstates to be degenerate. The set (1) decouples from the G — 1 states while the set (2) 
is annihilated by the G = states. In particular, the set (1) corresponds to the bound state 
wave-function found in the heavy quark limit of properly extended Skyrme type models M 



/ \ uiv) 

Wffc(Skyrme) = u {r) (f • r) al [r^Xh') = — {r ■ r {5 lh - r • cr lh ) X } a ■ 



(22) 



The interpretation of this situation is obvious: If the background soliton field contains only 
G = states, as it is the case in Skyrme type models, the only solution will be given by 
the set (1). However, if also G — 1 states are present an alternative solution will exist. 



^The single quark wave-functions, g^, carry the dimension (energy) 3 / 2 . 
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The soliton dynamics, from which the spectral functions fo and fx are computed, then 
determines which of these two modes, (1) or (2), leads to a larger binding and therefore has 
to be interpreted as the lightest baryon containing a heavy quark. 



IV. NUMERICAL RESULTS FOR THE HEAVY BARYON 

In analogy with the meson sector we have expanded the regularization function, Ra{uj, e), 
(|TT|) up to fifth order in the frequency uj. In figure 1 the resultant frequency dependence of 
the spectral functions ( ^0[) and (|2lD is shown for the special case m = 400MeV. As a result 
of the heavy quark transformation (0) the threshold for the (unphysical) decay of the heavy 
meson into a quark-antiquark pair is given by the energy of the lowest accessible light quark 
state. Therefore a singularity appears at u « — e va i for the mode coupling to the G = 
channel, while the radius of convergence is much larger in the case of the G = 1 one channel. 
The latter mode apparently leads to a smaller eigenfrequency, uq. 

The total energy of the system consisting of the soliton and the bound meson is given 
by Mq + E tot + ujq. As the total energy of the heavy meson in the trivial background 
is Mq + AM the binding energy of the heavy baryon becomes Eb = AM — uq. This 
binding energy is measured with regard to a decay into a soliton and a heavy meson. When 
increasing the heavy quark coupling constant, G$, the curves shown in figure 1 get shifted 
downwards by a constant amount. Hence the roots decrease and the binding energies of the 
baryon states containing a heavy quark become larger when the coupling gains strength. 
Obviously, a heavy baryon constructed from an eigenmode of type (2) is more strongly 
bound than a type (1) heavy meson. In table |lj the predicted roots, ujq, as well as the 
associated binding energies, Eb, of the heavy baryon are presented. For the modes of type 
(1) the eigenfrequency acquires a maximum for m ~ 500MeV. Since AM increases with 
the constituent quark mass, m, the binding energy for the associated heavy baryon exhibits 
almost no variation for m < 500MeV. When m is further increased the eigenfrequency 
decreases leading to a larger binding energy. This behavior of u is due to the decrease of 
e va i when m grows. On the contrary, the eigenfrequency corresponding to the eigenvectors 
of type (2) appears to be a quickly, monotonously decreasing function of m. This causes the 
binding energy of the corresponding heavy baryon to significantly increase with m. Since 
studies in the light quark sector of the NJL model favor 400MeV < m < 450MeV ITS 
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the results displayed in table [TJ suggest the interpretation that the NJL model predicts the 
baryon with one heavy quark to have a binding energy of about 560MeV in leading order 
of the 1/Mq expansion. A less strongly bound state is obtained possessing about half of 
that binding energy. As discussed above the latter state corresponds to the one which is 
identified as the most strongly bound baryon in Skyrme type models. In this context it has 
to be remarked that a state associated with the mode (2) does not appear in Skyrme type 
models because in those models the background soliton is a pure G = configuration. If the 
NJL model soliton had contained no G = 1 components the most strongly bound baryon 
would have been of the same structure as in purely mesonic soliton models. 

As in the bound state approach Jll| heavy baryon states with good spin and isospin are 



generated by canonically quantizing the collective isospin rotation, A(t), of the meson fields. 
This adds a term of the form Citr(xx' f A' ( (t)A(t)) t° the action. The index i = 1,2 refers 
to the distinct sets in eqs (pO| , pT|) . The time derivative of the collective rotation defines the 
angular velocity uj via A\t)A(t) = (i/2)u} ■ t. Empolying (iso) rotational invariance these 
coefficients can hence be computed from the matrix elements^ (/i|r 3 |z/) (v\£V W (3H£l\n) . As 
argued above H only couples to quark states with G — 0, 1. The bound state of set (1) has 
non-vanishing matrix elements only for quark levels with G = 0. In that case the matrix 
element of r 3 vanishes, hence C\ = 0, in agreement with the Skyrme model results ||. For 
the set (2) there is indeed a coupling when both quark levels and |z/)) are from the 
G* = 1+ channel. We find C 2 = C 2 [(C + 2)(C + 6)]/[(C + 2) 2 + 2], where the positive definite 
denominator is due to the normalization associated to the metric induced by the coefficient 
matrix of f\ in eq fll8|) . Although the collective rotation removes the degeneracy in ( we will 
argue that C 2 is negligibly small without going into the details of the calculation. For quark 
states and from the G* = 1 + channel the above mentioned matrix elements are 
similar to those entering the evaluation of the moment of inertia a 2 [22|. As the collective 
quantization gives rise to a factor 1/a 2 , C 2 can be estimated by the contribution of the 
= 1 + channel to the total moment of inertia. This contribution is very small because 
a 2 is dominated by the valence quark level, which resides in the G* = + channel, e.g. for 
m = 400MeV we find a contribution of only 2%. Hence we conclude that there is (almost) 



5 For the case of the bound state approach to hyperons these matrix elements have thoroughly 
been discussed by Weigel et al. in ref p 
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no hyperfine splitting in the heavy quark limit. For a finite heavy quark mass this conclusion 
should change because H will couple to channels with G > 2 as well. 



V. CONCLUSIONS 



We have analyzed fluctuations of mesons with a heavy quark in the background of the 
NJL model chiral soliton in the leading order of the heavy quark mass expansion. We have 
succeeded in showing that at this order the associated Bethe-Salpeter equation possesses 
the bound state solution known from Skyrme model studies. We have furthermore seen that 
this Bethe-Salpeter equation contains an additional solution which arises from the coupling 
of the heavy meson field to vacuum quark states. Also, the novel bound state appears in 
degenerate pairs. The actual computation reveals that in the NJL model the novel states are 
more strongly bound than those configurations known from the Skyrme model. Although 
this certainly represents an interesting feature in its own, future studies are required to find 
out whether or not this novel state is an artifact of applying the well-established 1/Mq 
expansion to the NJL model. In comparison a few aspects of the Skyrme model calculation 
for finite Mq PJTD|] are worthwhile to be mentioned. First, the restriction to the leading 
order of the 1/Mq expansion has been found to be a somewhat crude approximation to the 
realistic cases causing states, which are predicted to be degenerate at leading order, to acquire 
distinct binding energies when corrections of subleading order are taken into account. We 
assume that a similar mechanism will remove the degeneracy of the novel states in the NJL 
model. Second, table (1) of ref || shows that the transition from infinite to realistic values 
of Mq is indeed continuous. In particular the order of bound states remains unchanged. 
Hence one is inclined to conjecture that the novel states persist for finite Mq. Third, not 
all bound states, which are predicted within the large Mq limit are permitted for finite Mq 
because the limits Mq — > oo and r — > do not necessarily commute. This respresents the 
only mechanism which would prohibit the existence of the novel states in the case of finite 
Mq. These issues certainly raise interest for generalizing the NJL model studies to finite 
Mq. Such studies are considerably more involved because it is crucial to first perform the 
heavy quark transformation and to subsequently regularize the functional trace to have a 
consistent interpretation of the model. 

In practice the binding energies of heavy baryons are measured with respect to the decay 
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into a nucleon and the lightest meson containing the corresponding heavy quark. The exper- 
imentally observed binding energies are ~ 613±50MeV and ~ 625MeV for the A#(5641±50) 
and Ac(2285), respectively ||. It should be noted that these baryons are degenerate in the 
heavy quark limit. In reasonable agreement with these data our numerical analysis yields a 
binding energy of about 560MeV for the most strongly bound heavy baryon when the only 
free parameter of the model is adjusted to reproduce the properties of the light baryons. 
The less strongly bound heavy baryon with a predicted binding energy around 250MeV 
may eventually be associated with Ac (2625). This baryon is bound by about (~ 285MeV) 
against the decay into a nucleon and a .D-meson. We would like to emphasize that the main 
purpose of the present paper was to discuss the novel coupling scheme between the heavy 
meson and the soliton rather than to achieve a precise agreement with the experimentally 
observed binding energies. 
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APPENDIX 



In this appendix we will briefly provide the expansions we have been using for the regu- 
larization functions, which enter the Bethe-Salpeter equations for the heavy quark meson. 

In case no soliton is present these functions are obtained from the quark-loop in the 



self-energy, U(v ■ p) (cf. eq (47) of ref [pH ) 

d 4 k 



tr 



HU(v-p)H 



-iN, 



c 



tr 



j> + m) HH 



(23) 



(2tt) 4 ((k - p) 2 -m 2 ){v-p + ie) ' 
Again, v denotes the velocity of the heavy quark inside the meson while p labels the mo- 
mentum of the meson. We refer to the literature P!3|,Ii| on the treatment of this quantity 
in Euclidean space and regularization of the quark loop. When v ■ p is smaller than the 
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quark-ant iquark threshold, which in the heavy quark limit is idential to the light quark 
constituent mass, the self-energy has the Taylor expansion 



u(vp) = £n (n) (o)(v P )' 



(24) 



n=0 



where the superscript refers to the deriavtive with respect to v ■ p. In the proper-time 



regularization scheme of ref [H|, which actually defines the model, the first five coefficients 
of the expansion (2JI) are 



n(°>(o) 
n«(o) 
n^(o) 
n< 3 )(o) 
n^(o) 
n (5) (o) 



N c m 2 
16tt 2 ' 



57T 



4^) + 



Ay/TT 



m 



m 



exp I — — I — 27r erfc 



A^/tt 



m 



m 



exp r a? 



7r erfc ( — 
A 



N, 



c 



2T 0, 



4 



16vr 2 

N c _ exp 
16m7r 2 | 3 

N c 



m 
A 2 " 



7r erfc 



rn 



16m 2 7r 2 

iVc f /_8_ 
16m 3 7r 2 I \ 15 



4 

3 + 2A 



■^j I + Trorfo I — 



m 

1+ aJ 



exp 
+ 



A 

/ m 2 
\ A 2 

2A J 



H — erfc 
4 



exp ( 1 - T crfe [ — 



7T 
— ( 

4 



(25) 



We then obtain the binding energy, AM, of the heavy meson by determining the root from 
the trunctated expansion and including the contribution from the purely mesonic part of 
the action, A^l, (El) 



2G, 



^n (n) (0)(AM) n = 0. 



(26) 



r 3 n =0 

The decay constant, of the meson with the heavy quark is obtained by coupling external 



electro-weak sources [Q. This then relates ju to the field normalization, yll^O) 



G 3 JUW(0)M H 



(27) 



where M# refers to the mass of the heavy meson. Using the B-meson parameters fg = 
180MeV and Mb = 5.3GeV yields the numerical results for AM and G3 listed in table |[ 

As the procedure leading to the regularization function R\(u,e) has already been ex- 
plained in section 3, it sufficiecs to just list the expansion up to order uj 5 
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TABLES 

TABLE I. The proper-time cut-off, A, the heavy quark coupling constant, G3, the root of 
the Bethe-Salpeter equation AM, and the heavy meson binding energy, Em as functions of the 
constituent quark mass m. 



m (MeV) 


350 


400 


450 


500 


600 


A (MeV) 


641 


631 


633 


642 


672 


G 3 (l(T 5 MeV- 2 ) 


1.53 


1.68 


1.82 


1.94 


2.17 


AM (MeV) 


327 


353 


369 


401 


449 


E M (MeV) 


23 


47 


81 


99 


151 


TABLE II. The roots, 


luq, of the spectral functions (| 


(]) and (21) for the heavy 


meson modes 


which couple to the grand spin zero and one quark channels, respectively. Eb 


= AM- 


-u>o represents 


the associated binding ener 


gy of the heavy baryon. All data are in MeV. 






m 


350 


400 


450 


500 


600 


Set (1) ujq 


94 


102 


113 


114 


86 


Eb 


233 


251 


257 


287 


363 


Set (2) u>o 


-128 


-182 


-223 


-261 


-341 


Eb 


455 


535 


592 


662 


790 
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Figure Caption 



Fig. 1: The spectral functions and their roots for the eigenvectors defined in eqs (p0| 



and (pip . The constituent quark mass has been taken to be m = 400MeV. 
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